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In this paper, in the particular case of a concave flux function, we 
are interested in the long time behavior of the nonlinear process as- 
sociated in [Methodol. Comput. Appl. Probab. 2 (2000) 69-91] to the 
one-dimensional viscous scalar conservation law. We also consider 
the particle system obtained by replacing the cumulative distribu- 
tion function in the drift coefficient of this nonlinear process by the 
empirical cumulative distribution function. We first obtain a trajec- 
torial propagation of chaos estimate which strengthens the weak con- 
vergence result obtained in [8] without any convexity assumption on 
the flux function. Then Poincare inequalities are used to get explicit 
estimates concerning the long time behavior of both the nonlinear 
process and the particle system. 

Introduction. In this paper, we are interested in the viscous scalar con- 
servation law with C 1 flux function —A 

2 

(1) d t F t (x) = ^-dxxFtix) + d x (A(F t (x)), F (x) = H * m(x), 

where m is a probability measure on the real line and H (x) = l{ x >o} denotes 
the Heaviside function. As a consequence, H * m is the cumulative distribu- 
tion function of the probability measure m. Since A appears in this equa- 
tion through its derivative, we suppose without restriction that ^4(0) = 0. 
According to [8] , one may associate the following nonlinear process with the 
conservation law: 

(2 ) \x t = X Q + oB t -J^A'{H*P s (X s ))ds, 
[ Vt > 0, the law of X t is P t , 
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where (-Bj)t>o is a real Brownian motion independent from the initial random 
variable Xq with law m and a a positive constant. The process X is said 
to be nonlinear in the sense that the drift term of the SDE depends on the 
entire law Pt of Xf. More precisely, according to [8], this nonlinear stochastic 
differential equation admits a unique weak solution. Moreover, H * Pt(x) is 
the unique bounded weak solution of (1). For t > 0, by the Girsanov theorem, 
Pt admits a density pt with respect to the Lebesgue measure on the real line. 

We want to address the long time behavior of the nonlinear process solv- 
ing (2) by studying convergence of the density pt (see [2] and [3] for a simi- 
lar study in a different setting). Since the cumulative distribution function 
x — > H * P s (x) which appears in the drift coefficient is nondecreasing, con- 
vexity of A is a natural assumption in order to ensure ergodicity. Then the 
flux function —A in the conservation law (1) is concave. 

In the first section of the paper, after recalling results obtained in [8], 
we show that trajectorial uniqueness holds for (2) under convexity of A. 
Then we introduce a simulable system of n particles obtained by replacing 
in the drift coefficient the cumulative distribution function by its empirical 
version and the derivative A' by a suitable finite difference approximation. 
When A is convex, existence and trajectorial uniqueness hold for this system. 
Moreover, we prove a trajectorial estimation of propagation of chaos which 
strengthens the weak convergence result obtained in [8] . Unfortunately, be- 
cause the empirical cumulative distribution function is a step function and 
therefore not an increasing one, this estimation is not uniform in time. 

The second and main section deals with the long time behavior of both 
the nonlinear process and the particle system. We address the convergence 
of the density pt of Xt by first studying the convergence of the associated 
solution H *pt of (1) to the solution with the same expectation of the 

2 

stationary equation ^-d xx F 00 (x) + d x (A(F 00 (x)) =0 obtained by removing 
the time derivative in (1). For this result, no convexity hypothesis is made 
on A. Instead, one assumes A(u) < for u G (0, 1), A'(0) < 0, A(l) = and 
A'(l) > 0. In contrast, to prove exponential convergence of the density of 
the particle system uniform in the number n of particles, we suppose that 
the function A is uniformly convex. This hypothesis ensures the existence 
of an invariant distribution for the particle system. In [14], a necessary and 
sufficient condition on the drift sequence is established for existence of the 
invariant measure and convergence in total variation norm for the law of 
the particle system at time t to this measure. In the present paper, the key 
step to derive quantitative convergence to equilibrium consists in obtaining a 
Poincare inequality for the stationary density of the particle system uniform 
in n. This density has exponential-like tails and therefore does not satisfy a 
logarithmic Sobolev inequality. So the derivation of the Poincare inequality 
cannot rely on the curvature criterion, used, for instance, in [5, 6, 12] or [13] 
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for the granular media equation. Instead we make a direct estimation of the 
Poincare constant using the specific analytic form of the invariant density. 
To our knowledge, our study provides the first example of a particle system, 
for which a Poincare inequality but no logarithmic Sobolev inequality holds 
uniformly in the number n of particles. 

Assumption. Throughout the paper, we assume that A is a C 1 function 
on [0,1] s.t. ,4(0) = 0. 

1. Propagation of chaos. 

1.1. The nonlinear process. Let us first state existence and uniqueness 
for the nonlinear stochastic differential equation (2). 

Theorem 1.1. The nonlinear stochastic differential equation (2) admits 
a unique weak solution ({Xt,Pt))t>o- Fort>0, Pt admits a density pt with 
respect to the Lebesgue measure on BL The function (t,x) i— > iJ * Pt(x) is 
the unique bounded weak solution of the viscous scalar conservation law (1). 
Moreover, 

(3) Vt > X t - X is integrable and E(X t - X ) = -A(\)t. 

Last, if the function A is convex on [0,1], (2) admits a unique strong solu- 
tion. 

Proof. The first and third statements are consequences of Proposition 
1.2 and Theorem 2.1 of [8] [uniqueness follows from uniqueness for (1) and 
existence is obtained by a propagation of chaos result]. 

According to the Yamada-Watanabe theorem, to deduce the last state- 
ment, it is enough to check that when A is convex, then trajectorial unique- 
ness holds for the standard stochastic differential equation 

dX t = adB t - A'(H * Q t {X t ))dt 

where (Qt)t>o is the flow of time-marginals of a probability measure Q on 
C([0, +oo), M). Since for each t > the function x \— > A'(H * Qt(x)) is nonde- 
creasing, if (X t )t>o and (Y t )t>o both solve this standard SDE, then \X t — Y t \ 
is bounded by 

\X - Y \ + f siga(X s - Y S )(A'(H *Q S (Y S ))- A'(H *Q s (X s )))ds, 
Jo 

and then by \Xq — Yq\ which concludes the proof of trajectorial uniqueness. 

Existence of the density pt for t > follows from the boundedness of the 
drift coefficient and the Girsanov theorem. To prove (3), one first remarks 
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that by boundedness of the drift coefficient, for each t > 0, the random 
variable Xt — Xq is integrable and 

E(X t -Xo) = - f E(A'{H*P s (X s )))ds 
Jo 

= - f [ A'( f P s (dy))p s (dx)ds. 
Jo JR \J-oc J 

For s > 0, since by the Girsanov theorem P s does not weight points, 

J R A>(J X ^P s (dy)y s (dx) = [A(H*P s (x))£™ = A(l). Q 

Corollary 1.2. Assume that A is C 2 on [0, 1]. Then the function H * 
Pt(x) is C 1,2 on (0, +oo) x R and solves (1) in the classical sense on this 
domain. 

Proof. By the Girsanov theorem, for to > 0, the law Pt of Xt admits a 
density with respect to the Lebesgue measure on R. Hence (t, x) t— > H * Pt(x) 
is a continuous function on (0, +oo) x R with values in [0,1]. According to 
[11], Theorem 8.1, page 495, Remark 8.1, page 495 and Theorem 2.5, page 
18, there exists a function u with values in [0, 1], continuous on [0, +oo) x R 
and C 1,2 on (0, +oo) x R such that 

VxGR, u(0,x) =H*P ta (x), 

2 

V(t,x) € (0,+oo) x R, ^«(t,x) = y9„ii(M) 

By the uniqueness result for bounded weak solutions of this viscous scalar 
conservation law recalled in Theorem 1.1, Vt >to, H * Pt(x) = u(t — to,x). 
The conclusion follows since to is arbitrary. □ 

1.2. Study of the particle system. For n E N*, let (tJn(i))i<j< n be a se- 
quence of real numbers. In this section, we are interested in the ?i-dimensional 
stochastic differential equation 

(4) dXl' n = a dB\ - a n ^ 1 {X ^< X *.« } ^ dt, X l ' n = X l , l<i<n, 

where (i? l )j>i are independent standard Brownian motions independent 
from the sequence (Xg)j>i of initial random variables. 

In the next section devoted to the approximation of the nonlinear stochas- 
tic differential equation (2), we will choose a n (i) equal to the finite difference 
approximation n{A{i/n) — A{{i — l)/n)) of A'(-). For this particular choice, 
the nondecreasing assumption made in the following proposition is implied 
by convexity of A. 
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Proposition 1.3. Assume that the sequence (a n (i))i<i< n is nonde- 
creasing. Then the stochastic differential equation (4) has a unique strong so- 
lution. Let (Y^' 11 , . . . , Y™' n ) denote another solution starting from (Yq 1 , . . . , Y n ) 
and driven by the same Brownian motion (B , . . . , B n ). Then 

n n 

(5) a.s.yt > J2( X t n ~ Y r? < " Yof- 

i=l i=l 

In addition, if the initial conditions (Xq, . . . , Xq) and (Yq, . . ., Y n ) are s.t. 
a.s., Vi € {1, . . . , n}, X* < Yj (resp. X l < Y$), then 

(6) a.s., Vt > 0,Vt G {1, ... ,n} X\ n < Y/' n (resp. X^ n < Y/' n ). 

Existence of a weak solution to (4) is a consequence of the Girsanov the- 
orem. Therefore, according to the Yamada-Watanabe theorem, it is enough 
to prove (5) which implies trajectorial uniqueness to obtain existence of a 
unique strong solution. To do so, we will need the following lemma. 

Lemma 1.4. Let (a(i))i<i< n and (b(i))\<i< n denote two nondecreasing 
sequences of real numbers. Then for any permutation r € S n , 

n n 

(7) Ea(0&(r(0)<£o(.>(i). 

i=i i=i 

Proof. For n = 2, the result is an easy consequence of the inequality 

(o(2)-a(l))(6(2)-6(l))>0. 

For n > 2, we define t\ as r if r(l) = 1 and as r composed with the trans- 
position between 1 and t _1 (1) otherwise. This way, t\(1) = 1. In addition, 
using the result for n = 2, we get Y^=i a (*)K r (*)) — Ya=i a(i)b(Ti(i)) . 

For 2 < j < n — 1, we define inductively Tj as t 3 -_i if t 3 -_i(j) = j and as 
Tj_i composed with the transposition between j and rf_^(j) otherwise. This 
way, for 1 < i < j, rj(i) = i. Again by the result for n = 2, one has 

TL Th Tb Th 

J2 a(i)b(r(i)) < £ a(i)b(n(i)) < £ «(0&fa«) < • ■ • < £ a(i)b(r n ^(i)). 

i=l i=l i=l i=l 

We conclude by remarking that r n _i is the identity. □ 

We are now ready to complete the proof of Proposition 1.3. 

Proof of Proposition 1.3. Let (X 1 ' 11 , . . .,X n ' n ) and (Y 1,n , . . . , Y n ' n ) 
denote two solutions. The difference 

n n 

j2(xr-Yrf-j2(xi-Ytf 

i=i i=i 
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is equal to 

( 8 ) 2 J Q ~ Y s' n )(^n\^l{ Y J- n <Y*< n }j ~ a n(^2\xi' n <xi' n }j^j ds ' 

By the Girsanov theorem, for any s > the distributions of (Xg ,n , . . . , X™' n ) 
and (Yg' n , . . . , YJ 1 ' 71 ) admit densities w.r.t. the Lebesgue measure on M. n and 
therefore dF®ds a.e. the positions X]> n , . . . , X™< n (resp. Y s 1 ' n , . . . , YJ 1 '™) are 
distinct and there is a unique permutation G 5 n (resp. r s y € <S n ) such 

that X s s v >' <X s sK1 <---<X s s() (resp. Y s " K <Y S " K <■■■< 

Y rJ{n),ny There f ore (Tp ^ a e . ; 

- {Y l{Y-j> n <Y? n }j ~ ° n 1 {X J s ' n <Xl n }j ^ 

is equal to 

n 

J2an(^)((xf^ n - YJ?®' n ) - (Xf^' n - Y^^ n )). 

i=l 

The sequence (a n (i))i<i< n is nondecr easing. Applying Lemma 1.4 with 6(i) = 

xf {{) ' n and r = (rf [- 1 o T J then with = Y? r ®> n and r = (rf J" 1 or* , 
one obtains that the integrand in (8) is nonpositive cflP (8) ds a.e. Hence (5) 
holds. 

Let us now suppose that a.s. Vi G {1, . . . ,n}, Xq < Yq and define v = 

inf{t > : 3i G {1, . . . , n}, X^' n > Y^' n } with the convention inf = +oo. From 

now on, we restrict ourselves to the event {v < +oo}. Let i G {1,. .. , n} be 

such that Y"J' n = X*' n . There is an increasing sequence (sfc)fc>i of 

positive times with limit v such that V/c > 1, Qn(X)?=i l/yj'. n <-v*> n i) < 

^- s fc — s /c j 

). Since (a n (i))i<j< n is nondecreasing, by extracting a 

subsequence still denoted by (sfc)& for simplicity, one deduces the existence 
of j G {1, . . . ,n} with j ^ i such that VA; > 1,X*£ < and Y^> < y^ n . 
Since s k < v, X l s £ < X j s '™ < Y^ n < Y^ n . By continuity of the paths, one 
obtains X]] n = Xl ,n = Yj' n = Y*' n . Now since the probability of the event 



3* l5 i 2 , i 3 dist. in { 1 , . . . , n} , 3t > X l L + aB] 1 = X l 2 + aB\ 2 = X l 3 + aBf 

is equal to 0, the Girsanov theorem implies that a.s. V/ G {1, . . . ,n}\ {i, j}, 
X l u > n £ Xl> n = Xy n . In the same way, Y^ n / Y^ n = Yj> n . By continuity of 
the paths and definition of v one deduces that for k large enough, and for 
every t G [s k ,u], 

n n n n 

Y l {Y t l ' n <Y t l ' n } - Y 1 {X l t ' n <X l t ' n y Y 1 {Y^' n <Y t J ' n } - Y l {X l t - n <x{' n }- 
1=1 1=1 1=1 1=1 

l&j l ^i,3 l&j tyhj 
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Since a.s. dt a.e., Y^' n / Y/' n and (a n (i))i<j< n is nondecreasing, one obtains 
that a.s. dt a.e. on [sj.,i>], 

+ a n ^2 1 {Y t l - n <Y t j - n }j 

< a n 1 {x l t ' n <x{- n }j + Qn (Y1 1 {x l t - n <x l / 1 }j ■ 

By integration with respect to t on [s^, v], this implies that a.s. Y^ n — X]) n + 
Yj' n - Xl' n > Y l s - k n - X l s £ + y/> - A7> > 0. Therefore ¥(v < +oo) = 0. 

When a.s. for i G {l,...,n}, Xq < Yq, one obtains that for e > the 
solution (y n,£ , . . . , Y t n ' n ' £ ) to (4) starting from (Yq 1 + e, . . . , Y n + e) is such 
that 

a.s., Vt > Vi G {1,. .. ,n} X*' n < y/' n ' £ . 
Since by (5), Y t l,n ' e < Y t l,n + y^e, one easily concludes by letting e — > 0. □ 

1.3. Trajectorial propagation of chaos. From now on, we set 

(9) VnGN*,ViG{l,...,n} a n (i) = n( A[ - A 



n l \ n 



and assume that the initial positions (Xq)j>i of the particles are independent 
and identically distributed according to m. We prefer to define a n (i) with 
the above finite difference approximation of the choice A'{ijn) made in [8] 
because the sum Y^i=i a n('i) which plays a role in the long time behavior of 
the particle system is then simply equal to nA(l). One could also obtain 
trajectorial propagation of chaos estimates similar to Theorem 1.5 below for 
the choice a n (i) = A'(i/n). 

In the present section, we also suppose that A is a convex function on 
[0,1]. By Theorem 1.1, for each i > 1, the nonlinear stochastic differential 
equation 

(10) | A = A + aB\ - jT A'(H * P S (XD) ds, 

lvt>0,the law of Xj is P t , 

has a unique solution and for all t > 0, the law Pt of X\ does not depend on 
i. Under a Lipschitz regularity assumption on A' , we obtain the following 
trajectorial propagation of chaos estimation. 

Theorem 1.5. // A: [0, 1] — ► M is convex and A' is Lipschitz continuous 
with constant K, then 

kH 2 



Vn> 1,V1 < i < n,Vt > E( sup (X l s ' n - X l s ) 2 ) < 

se[o,t] 



(yn 
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Proof. Let us write Ei=i( x t - x t) as 



L {xi<xi} 



r t n ( ( n 

:/ o £w»-*i)(.„(x;i 

ft n 

+ 2/ Y.( X T- X l)C{s,Xl...,X™)ds 
Jo ~1 



{Xl> n <Xl> n } 



ds 



a(±Ti 1 
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where C(s, . . . , X") is equal to 

A\H*P s (Xl))-n^±l {x ^ Xls} 

Like in the proof of trajectorial uniqueness for (4), because of the convex- 
ity of A, the first term of the r.h.s. is nonpositive. Moreover, by Lipschitz 
continuity of A' , 

\ 2 



n t- 



A (^ n EVj<*j} n 



/ 1 n — 1 



n +(f*p s (x:)-0) d<9. 



For s > 0, as the variables X\ are i.i.d. with common law P s which does not 
weight points and H * P s (X l s ) is uniformly distributed on [0, 1], 

2\ 



jT 1 E^^(fr*p a (xi)-i {X i< xi} ) + (fr*p a (^)-0) 



de 



£ H(H * Ps(Xl) - l {x| < xn ) 2 ) + /' E((F * Ps(Xl) - Of) dO 

■ i ■ J 



= (n - 1)E((H * P s (Xi))(l - * P s (Xi))) + 1/6 
= ra/6. 

Using the Cauchy-Schwarz inequality, one obtains 

rt 



E sup £(*j«-*j) 2 <2 
V«e[o,t] i=1 / 



\ 6ra 



(is 



< 



2K 



E sup £(Xi' n -X*) 2 da. 



VE JO \ \u€[0,s] i=1 
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By comparison with the ordinary differential equation a'(t) = 2K\J —jp-, one 
concludes that 

( n \ K 2 f 2 

Vi>0 E sup -*j) 2 <-— . 

Exchangeability of the couples ((X l ' n ,X l )) ie r x n \ completes the proof. □ 

Remark 1.6. One could think that assuming that A is uniformly con- 
vex: 

(11) 3a>0,V0<x<y< 1 A' (y) - A' (x) > a(y - x) 

would lead to a better estimation. Indeed, then for every i G {1, . . . ,n — 1}, 

Ai+l)/n 
li/n 



a n {i + l) -a n (i) = n 



A'(z)-A'(z-^ 



a 

ax > — . 

n 



But since even in this situation, the nonpositive term 



^2( x l' n x l)y a ^[^2 1 {xi<xi}J a n\Yl l {xi' n <xi' n }J j 

vanishes as soon as the order between the coordinates of (Xl' n , . . . ,X™ ,n ) is 
the same as the order between the coordinates of (Xj, . . . , X™), we were not 
able so far to improve the estimation. 

Corollary 1.7. Under the hypotheses of Theorem 1.5, let m be a 
probability measure on R such that Vx G M, H * rh(x) < H * m(x). If for 
some random variable U\ uniform on [0,1] independent from (B l )i>i, Xq = 
ml{x : H * m(x) > U±} and (Y^ 1 )t>o denotes the solution of the nonlinear 
stochastic differential equation 

(12) | Yt 1 = Yi + oB\ - jf* A'(H * P^Ys 1 )) ds, 

{ Vt > 0, the law of is P t , 
with Yq = inf{x : H * rh(x) > U±}, then 

p(vt>o,x i 1 <y i 1 ) = i. 

Moreover Vt > 0, Vx G R, H * Pt(x) < H * Pt(x). Last, the function 1 1— > 
EjYj 1 — X}\ is constant. 

Remark 1.8. At least when m and rh do not weight points, one has 
a.s. A'(H * Po(X^)) = A'(H * P (Y£)) since H * m(X^) = H * fh(Yj) = Ux. 
Therefore a.s. d{Y x -X 1 )q = Q and one may wonder whether a.s. Y^ — 
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X} does not depend on t. If this property holds, necessarily, a.s. dt a.e. 
A'(H * P t (X})) = A'(H * PtiYt 1 )). If A' is increasing, a.s. for alH > 0, H * 
Pt(Xj) = H *p t (Y t 1 ) with p t and pt denoting the respective densities of Pt 
and P t . If A is C 2 , the Brownian contribution in d(H *p t (X^) — H * p t (Y^)) 
given by Ito's formula vanishes, that is, pt(X}) = p t (Y t ) and Vu G]0, 1[, 
p t ((H*p t )- 1 (u)) =p t ((H*p t )- 1 (u)) or equivalents ((H *p t )~ 1 )'{u) = ((H * 
pt)~ l )'(u). Hence Yj 1 = X\ +c for a deterministic constant c which does not 
depend on t according to (3). Letting t — > 0, one obtains Yq = Xq + c. This 
necessary condition turns out to be sufficient as {X} + c)t>o obviously solves 
the nonlinear stochastic differential equation (2) starting from Xq + c. 

Proof of Corollary 1.7. For (C/j)j> 2 a sequence of independent uni- 
form random variables independent from (U\, (i? l ),>i), we set 

Vi>2 X l = inf{x : H * m(x) > U t } and Yq = inf{x : H * rh(x) > Ui}. 

Since H*rh < H*m, a.s. Vi >1,Yq> Xq. From Proposition 1.3, one deduces 
that the solutions (X t ' n , . . . , and (Y t ' n , . . . , Y[ l,n ) to (4) respectively 

starting from (Xq , . . . , Xq) and (Yq, . . . , Y n ) are such that 

a.s., Vn>l,Vt€{l J ... ) n},Vt>0 Yj' n > X l t ' n . 

Since, by Theorem 1.5, for fixed t > 0, one may extract from (X t ' n , Y t ' n ) n >i 
a subsequence almost surely converging to (X^Y^), one easily deduces that 
P(Vt > 0, X t x < Y^) = I. Hence 

Vt > 0, Vx G R H* P t (x) = P(y 4 1 < x) < Fix} <x) = H*P t (x). 

Since |Y/ - - l^ 1 - = Y^ 1 - Yq 1 - (X} - X^), (3) ensures that E\Yf - 
X}\ S [0, +oo] does not depend on t. □ 

2. Long time behavior. In this section we are interested in the long time 
behavior of both the nonlinear process and the particle system. According 
to (3) and the equality J2i=i a n(i) =nA(l) which follows from (9), we have 
to suppose A(l) = in order to obtain convergence of the densities as t 
tends to infinity. We address the convergence of the density pt of X t by first 
studying the convergence of the associated cumulative distribution function 
Ft under the following hypothesis denoted by (H) in the sequel: 

A(0) = A{1) = 0, A'(0)<0, 

(H) 

4'(1)>0 and VuG (0,1) A{u) < 0. 

These assumptions determine the spatial behavior at infinity of the drift 
coefficient in (2). 

To prove exponential convergence of the density of the particle system 
uniform in the number n of particles, we make the stronger assumption of 
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uniform convexity on A. The key step in the proof is to obtain a Poincare in- 
equality uniform in n for the stationary density of the particle system. This 
density has exponential-like tails and therefore does not satisfy a logarith- 
mic Sobolev inequality. So the derivation of the Poincare inequality cannot 
rely on the curvature criterion, used, for instance, by Malrieu [12, 13] when 
dealing with the granular media equation. Instead, we take advantage of the 
following nice feature: up to reordering of the coordinates, the stationary 
density is the density of the image by a linear transformation of a vector 
of independent exponential variables. And it turns out that the control of 
the constant in the n-dimensional Poincare inequality relies on the Hardy 
inequality stated in Lemma 2.18 which is a one-dimensional Poincare-like 
inequality. To our knowledge, our study provides the first example of a par- 
ticle system, for which a Poincare inequality but no logarithmic Sobolev 
inequality holds uniformly in the number n of particles. 

2.1. The nonlinear process. In this section, we are first going to obtain 
necessary and sufficient conditions on the function A ensuring existence 
for the stationary Fokker-Planck equation obtained by removing the time- 
derivative in the nonlinear Fokker-Planck equation 



satisfied by the density of the solution of (2). Under a slightly stronger 
condition, the solutions satisfy a Poincare inequality. 

Lemma 2.1. A necessary and sufficient condition for the existence of a 
probability measure fi solving the stationary Fokker-Planck equation 



in the distribution sense is A(l) = and A(u) < for all u S (0, 1). Under 
that condition, all the solutions are the translations of a probability measure 
with a C 1 density f which satisfies 



(13) 





-^A(H*f(x)) and 



(14) 




If A' 



(0) < and A'(l) > 0, then 




when x 



when x — > +oo, 



— oo 
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(15) 



-a 2 



x dy I 2A'(0)/(x)' 

o f(y) ~ I ^ 2 

2A'(l)/(x)' 



iu/ien x — ► — oo, 
iu/ien x — ► +oo, 



and aZZ i/ie solutions satisfy a Poincare inequality and have a finite expecta- 
tion. Last, if the function A is C 2 on [0, 1], then f is C 2 and satisfies 

(16) f"{x) = ~- 2 A"(H * f(x))f 2 (x) + Q^. 

Proof. Let /U be a probability measure on R solving the stationary 
Fokker-Planck equation. The equality \d xx [i = —d x (A'(H*fj,(x))n) ensures 
that \x does not weight points. Hence the stationary equation is equivalent 
to d xx { 2 2^i + A(H * fJ>(x))) = 0. One deduces that [i possesses a C 1 density 
/ such that 

(17) VxGR f(x) = --^A(H*f(x)) + ax + (3, 

a z 

for some constants a and /3. Since A(0) = 0, letting x — ► — oo then x — > 
+oo in the last equality, one obtains a = (3 = A(l) = 0. For u G (0, 1), since 
u = H * f(x) for some x G R and H * f is not constant and equal to u, 
the Cauchy-Lipschitz theorem and (17) imply that A{u) / 0. Since / is 
nonnegative, A(u) < 0. Hence A(l) = and A{u) < for all u G (0, 1) is a 
necessary condition. 

Under that condition, a probability measure /i solves the stationary Fokker- 
Planck equation if and only if its cumulative distribution function H * /x(x) 
is a C 2 solution to the differential equation 

(18) (f/(x) = —zA(<p(x)), xGM. 

By the Cauchy-Lipschitz theorem, for each v G [0, 1] this equation admits 
a unique solution tp v defined on R with values in [0, 1] such that (p v (0) = v. 
Moreover, as A(0) = A(l) = 0, (po = and ipi = l and 

(19) Vv£ (0,l),VxGR 0<<p v (x)<l. 

For v G (0, 1), since cp v is nondecreasing and ip v (x) = v — \ Jq A(ip v (y)) dy, 
necessarily lim y ^ +00 tp v (y) = 1. In the same way, lim^-oo <p v (y) = and ip v 
is an increasing function from R to (0, 1) with inverse denoted by c^" 1 . The 
uniqueness result for (18) implies that Vi> G (0, 1), Vx G R, ip v {x) = <P\/2{x + 
92^2 (u)). Therefore the solutions to the stationary Fokker-Planck equation 
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are the probability measures obtained by spatial translation of the proba- 
bility measure with density f{x) = ^^(x) which satisfies (14) according to 
(18). 

Let us now suppose that A'(0) < and A'(l) > 0. When x — > +oo, 
f(x) = --A\l-j^ f(y)dy)~-£±J x f(y)dy. 

B y ( 14 )' 7$ = Q°gf( x )y = -3t A '(<Pi/2( x )) converges to -^P- as x -> 
+oo. This implies that los ^J x ^ converges to — 2A JP and that xf{x)lt x >Q\ is 
integrable. Moreover, since J +o ° = +oo, ~ ^TJ Jo "fS dy ~ 

2 

2A'(i)/(x) ' as 00 ~ > +°°- I n the same way, one obtains the equivalents given in 
(15) when x — ► — oo and checks the integrability of the function x/(x)l| a ,< }. 
From (15), one has 



lim I f{y)dy f 

00 J — oo a; 



dy a 4 



/(y) 4(A'(0))2 



.r- 

and 

rx dy a 4 



r+oo r 

lim / f{y)dy \ 



2 ' 



o f(y) 4(A'(1)) 

By Theorem 6.2.2, page 99 of [1], one concludes that the measure with 
density / satisfies a Poincare inequality. 

By (14), the function / is C 2 as soon as the function A is C 2 on [0,1]. 
Moreover, f"(x) = -^A" (H * f(x))f 2 (x) - ^A'(H * f(x))f (x) which com- 
bined with (14) implies (16). □ 

Remark 2.2. When A is a C 1 convex function on [0, 1] such that A(0) = 
A(l) = and A'(u) < for some u € (0, 1), then the necessary and sufficient 
condition in Lemma 2.1 is obviously satisfied. Since (14) implies 

(log f(x))" -( f 'W-( - 2 A^'(tf * /(*))/(*) \ ' 



f(x) J V f{x) 
= ~A"(H*f(x))f(x)<0, 

the probability measures solving the stationary Fokker-Planck equation 
admit log-concave densities with respect to the Lebesgue measure. Log- 
concavity is a property stronger than the existence of a Poincare inequality 
(see [7]). 

Example 2.3. Using (18) and (19), the following two choices for A lead 
to exact computations and different tails for the stationary densities: 
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• if A(x) = \x{x - 1), one gets log( 1^7^) ) = ^A^' that is ' 

^ 2(X) = IT^ ^ ^ 2(x) = 4^cosh'(,/2. 2 ) ; 

• if A(x) = x 3 — x = x(x — l)(x + 1), 

1 2e~ ix / u2 



When ^4(1) = and A(u) < for all (0,1), a natural question is how 
to link the translation parameter of the candidate long time limit of the 
marginal Pt solving the stationary Fokker-Planck equation to the initial 
marginal m. When J R \x\m(dx) < +00, by (3), for all t > 0, E(A' t 1 ) = E(Xq). 
Therefore the translation parameter is chosen in order to ensure that the 
invariant measure has the same mean as the initial measure m. 

Let us denote by pt the density of Pt and by Ft = H * Pt its cumulative 
distribution function. 

Theorem 2.4. Let A be C 2 on [0,1] satisfying (H). Assume that m 
admits a density po such that J R \x\p (x) dx < +00 and J R {po{x) ~J{x) X)? dx 
is small enough where p^ denotes the stationary distribution with same 
expectation as pq. Last, we suppose that A andpo are such that p is a smooth 

solution of (13). Then f R feiMgg^ dx converges to exponentially fast 
as t — > +00. 

By a smooth solution of (13), we mean that p possesses enough regularity 
and integrability so that the formal computations made in the proof below 
are justified. 

Example 2.5. When A(x) = \{x 2 — x), one easily checks that the func- 
tion 4>(t,x) = —F t (x + |) solves Burgers' equation 

dt<t> = ^rd xx <j) - ^d x <p 2 , 0(0, x) = -F (x). 

By the Cole-Hopf transformation, ip(t,x) = exp(— ^ Pf^ (j)(t, y) dy) solves 
the heat equation 

d t Tp = Y d xx'>P, ijj(0,x)=exp(J^ J F (y)dyj. 

Since F t (x) = <r 2 %^(t,x — |), one deduces that 



(20) F t (x) 



k er^^^Foiy)^, y) dy/jaV^t) 

J R e -(x-t/2-y)y2vH^ y ) dy /(ay/tort) 
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If x denotes the expectation associated with the cumulative distribution 
function Fq, one has Jf^ Fq(z) dz = f^~°°(l — Fo(z)) dz. Since 

/X fX fX 

F (z)dz= / F (z)dz- / (l-F (z))dz+(x-x), 
-oo <j — oo Jx 

one deduces that the function $(0, x) = e~^ x ~ x ^ a3 i/}(0, x) [resp. ip(0, x)) is 
bounded on W + (resp. R_) and converges to 1 as x tends to +00 (resp. —00). 

Let us deduce the limit of Ft(x) as t -> +00. Writing the integral for 
y S R as the sum of the integrals for y S M_ and for y £ R + , and making the 

change of variables z = V ~*^ 2 (resp. z = V ~^~/j^ 2 ) i n the first (resp. second) 
integral, one obtains 



Jr 



-z 2 /2 



l {z<Vt/(2a)-x/(aVt)} 



aVtz + x - ^jip(o, aViz + x - —^= 



+ e (--)/- 2 / e -#/*i 



{z>-Vt/{2a)-x/(aVi)} 



x Fq ( a Viz + x + -^\ij}(o, a\ftz + x + ^ 



'2tt 

By the Lebesgue theorem, the first term of the right-hand side converges to 

whereas the second term converges to ei x - x )/ a . Replacing Fq by 1 in the 
above computation, one obtains that the denominator in (20) converges to 

1 + e (*-20/<x 2 _ Therefore 

p (x-x)/a 2 

VxGE lim F t (x) 



(AToo' 1 ^'' 1 + e {x-x)/a 2 ' 

Notice that in the same way, one may also obtain the limit of the density 
J R ((y + t/2 - x)/(aH))e-(*- t / 2 -^ 2 /^F (y)m y) dy/ja^t) 

Pt{X) j u e-^-y^ 2 Mo,y)dy/(aV^i) 

1 (f R e-^- t / 2 -^ 2 /^F (y)my)dy/(crV2rt)\ 2 



a 2 



( k e-W-vrw*)F (yM0, y) dy/(a^t) \ * 
V J m e-^ t / 2 ~y) 2 /^ 2t mO,y)dy(aV2^i) J ' 
One easily checks 

^ / e (x-x)/a 2 e 2(x-x)/a 2 

VxGR lim pt(x) 



t-^+oo^ > a 2 y 1 + e (x-x)/a 2 (l + e (x-x)/a 2 y 

1 



4a 2 cosh 2 ((x - x)/2a 2 ) 
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In order to prove Theorem 2.4, we are first going to check exponential 
convergence of F t to the cumulative distribution function of p^. Let 
Gt = F t — -Fqo . Since for a random variable X with cumulative distribution 
function F, K(X) = Jd~°°(l — F(x))dx — J^_ 00 F(x) dx, the equality of the 
expectations associated to Ft and writes J R Gt(x) dx = 0. This very con- 
venient expression of the link between pt and poo is one main reason for first 
considering the convergence of Gt to 0. In order to prove this convergence, 
we need the following result. 

Lemma 2.6. Under the assumptions of Theorem 2.4, one has 

, \ f G 2 (x) , r fGAx)\' 2 . . , 

(21) / l-^r) Poo(x)dx 

where c denotes the constant in the Poincare inequality satisfied by p^. 
Moreover 

(p t (x) -Poo{x)) 2 



(22) 



and 



Poo(x) 
G t (x) 



■ dx 



Poo{x) 



Poo {x)dx + \ I G t {x) 2 A"(F 00 )(x)dx 
cr z Jm. 



(23) [^ldx<c[ Mx) -fj x))2 dx. 

JR Poo(X) JR Poo[X) 

Remark 2.7. When A is convex, (23) is a consequence of (22) and (21). 

Proof of Lemma 2.6. As J R G t (x) dx = 0, (21) is the Poincare inequal- 
ity satisfied by p^ written for the function Gt/p Q 



-'oo • 



Since f G 'M V- G tW Gt(x)p> (x) , 
hlI1Ce ^P^x)> ~ P^M ~ Poo (x)* ' OIle haS 



Poo(x)J JR Poc(x) JR Plo{x) 

Gl^p'^xf 



dx 



(Pt(x) - Poo (x)) 2 dx r G 2 t (x)p'^(x) dx 



Poc(x) JR Ploix) 

G 2 t {x)^{xf 



vlo{x) 

Since p^ solves (16), one easily deduces (22). 



dx. 
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Writing G\ (y) as 

/ ry r+oo 

2 [ 1 {y<o} Gt(pt ~ Poa)(x) dx - l{ y >o} / Gt(pt-Poo)(x)dx 

\ J — oo Jy 

one obtains 

(24) / < ^{x)dx = -2 f G t {pt-Poo){x) f —^dydx. 



Poo Jm. Jo PcoVV) 

By (15), and since is bounded from below and above on each compact 
subset of the real line, 



3C>0,Vxe 



dy 



o Poo(y) 



C 

< 



Poo(x) 



Using the Cauchy-Schwarz inequality in (24), and inserting the latter bound, 
one obtains 



Poo \JRPoo J \JR Poo{x) 

One easily deduces (23). □ 

According to (23), the exponential convergence of J m {Mx) ~J^ {x))2 dx to 
zero is a stronger result than the exponential convergence stated in the next 
lemma. 

Lemma 2.8. Under the assumptions of Theorem 2.4, there is a positive 
constant C such that if J K —^(x) dx is small enough, then 

r C 2 p~ ct r r 2 

Vt>0 / ^{x)dx<-— / ^(x)dx. 

JR Poo & JM. Poo 

Proof. According to (14), one has ^F^ + (A(F 00 )) / = which also 

writes ^ = —■^ I A'{F oa ). Combining these equations with (1), then using 
Young's inequality, one easily obtains for e > 0, 

1 d r G h ^ 

-— / —(x)dx 
2dtJ RPo y 



'oo 

■2 



o- 2 f (G t {x 



(25) 



Poo (x) dx 



2 Jk \Poo(x) 

f (A(F t ) - A(F 00 ) - J 4'(F 00 )G t )(x) ^ 
Jr 



Poo{x) 
ll\\2 r 



dx 



< 



( a 2 \ f (G t . . V , . , \A"f^ f Gf(x) , 
\ £ -^) X Poo(x)dx + ^-f^ / -^4-dx. 

V 2 / Jr \Poo J 16e JrPoo(x) 
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Since 



(26) 



\JR VPoo{X) V / 

Vr Poo (a;) 



|Gt| is bounded by 1 and Poo^'^-Foo) = ~-giA x -^■"(-^oo) is bounded, one 
deduces from (22) that 



\\G t \t <\\\AA"U I ^(x)dx+(lA f (^(x))' 2 poo (x)dx 

0" 4 JR Poo \ JR \Poo J 

Inserting this bound in (25) and using Young's inequality, one deduces that 
for rj > 0, 



1 d r G i, ^ 

— — / — (x)ax 



+ \\AA"\U\A"\\lf f G 



4ecr 4 Vr poo 

/2 X 2 II Ainu 



I — -(x) dx 

JR Poo 



-{ e+r] -Y)L(t {x) ) pco[x)dx 

, ( \\AA»\U\A»\\l \\A"\\l \ dx \ 2 



One easily concludes with (21) and Lemma 2.10 below. □ 

Remark 2.9. (i) After reading this proof, one may wonder whether one 
could replace the upper bound in (25) by 



-*L) [ (<k ( -^*- , \\ A "\\ 2 co f g 2 

2 / Jr VPoo 



V \\A" r r G z 

-(x) Poo { x )dx+^-^ / -L(x)dx 
J 16e Jr Poo 



using \\Gt \\oo < 1- If the constant c in the Poincare inequality (21) was smaller 
than r^p-, one could deduce exponential convergence of L— L (x)dx to 

II II oo /^oo 

even for large values of J R ^-(x) dx. In case A(x) = ^{x — x) (see Example 
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2.5), one has H-A"!^ = 1 and 

>4 C t 4 / y 2 e- 2 y<ly = a 4 



\A»\\ 2 ' 
J^ 1 II OO 

and this approach does not work. 

(ii) Convexity of A implies nonnegativity of the term A(Ft) — A(F OQ ) — 
A'lF^Gt which appears in the right-hand side of the first displayed equality 
in the proof. One may wonder if one could exploit this property to obtain 
exponential convergence of p t to p^ even if po is not close to poo. We have 
not been able to do so. 

Proof of Theorem 2.4. By (14), p'^ = —^A^F^p^ and \\poo\\oa < 
2 ^|°° . The Fokker-Planck equation (13) for pt ensures that 

1 d f (pt(x) -Poo(x)) 2 dx 



2dtJu Poo{x) 



4 / (f w) 

2 Jm. \Poo J 

(A!{F t )-A\F 00 )){x)( Pt -p 00 ){x)(^{x) S \ dx 

\Poo J 

(A' (F t ) - A' (F 00 ))(x)p 00 (x) (—(x)\dx. 

\Poo J 



Then, using Young's inequality and (26), one easily checks that for e,rj > 0, 
Id f (pt(x) -Poo(x)) 2 _ dx 



2dtjR Poo(x) 

- ( v+£ ~ t) X (p^^O p °°^ dx 

4e7R Poo X 



+ / (A / (F t )(x)-A / (F 00 )(x)) 2 p 00 (x)dx 
477 Jr 

*(* +e -TU(£ (a °) Poo(:c)dx 

4e \ii Poo(^) 
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+ \\A^ x Aj\A^ r Gl {x)dx ^ 

4n (J 4 Jr Poo 

By (23) and Lemma 2.8, for / R (po( ^~^ (x))2 dx small enough, the last term 
of the r.h.s. is smaller than / R (po( ^~^ (x))2 dx. Since J R (^(x))' 2 P oo(x) dx 

O-Poo(z )) 2 J 

Poo (a) 



is greater than i J K (p*^) dx, one easily concludes by Lemma 2.10 be 



low. □ 



Lemma 2.10. Assume that u : R+ — ► M+ satisfies 

Vt > ^(i) < /3n(t) («(t) - a) + 7e" 5 * 

/or some constants a,/3, S > and 7 > 0. 
7/7 = and n(0) < a, iden 

an(0)e" Q/3t 

Vt > n(t) < 



a + n(0)(e" a ^ - 1)' 

// n(0) < f and 7 < £den u(t) converges to exponentially fast as 
t — > +00. 

Proof. When 7 = 0, as long as it(t) € (0, a), one has 

du, J 1 1 \ 

dt \u(t) a — u(t) J 

and after integration one obtains the desired estimation. Since the upper 
bound is not greater than n(0) and u(t) = => Vs > t,u(s) = one easily 
concludes. 

Now when 7 G (0, ^j-), one has /3a (a — a) = 7 for some a G (0, f ) and 

d / . a \ + ^ du . . 

— I «(t) A — — a I = l{ a <u(t)<a/2} -^(*) < 0. 

Hence when u(0) < § , \/t > 0, n(t) < n(0) V a and 

^(t) < -P(a - tt(0) V a)u(t) + 7e~ 5 *. 
For = e /3("-«(0)va)t n( -^ Qne Educes 

^( t ) < 7e (/8(o-u(0)Vo)-«)t 

and one concludes by integration of this inequality that n(t) is bounded by 

C(l + t) e -[(/3(«-«(o)Vo))A5]t i n 
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2.2. The particle system (4). Let us suppose that A(l) = and that 
the first-order moment associated with the initial probability measure m 
is defined and equal to x. As in the case of the granular media equation 
considered by Malrieu [12, 13], the direction (1, 1, ... , 1) is quite singular for 
the particle system. Indeed, 

n 

d{xi' n +...+xr)=°Y. dB i 

i=l 

which prevents the law of (Xl' n , . . . , X™ ,n ) from converging as t — > +oo. 
Following [12, 13], one introduces the hyperplane M. n = {y = (yi, . . . ,y n ) G 
K n : yi + • — VVn = nx} orthogonal to this singular direction and denotes by 
P the orthogonal projection on M n and by P the orthogonal projection 
on {i/ = (i/ 1 ,...,!/ n )£M":t/i + - + ^ = 0}. Since E"=i a n (i) = n(A(l) - 
A(0)) = 0, the orthogonal projection (Y/' n = x + X t i,n - ± £™ =1 ^/' n )i<i<n 
of the original particle system on is a diffusion on this hyperplane solving 

(27) dy/' n = a^—± dflj - - £ dB? - a n ( £ 1 { y/-<y t - } ) 

Propagation of chaos for the projected system is a consequence of the 
following estimate. 

Proposition 2.11. Assume that A is convex, such that A' is Lipschitz 
continuous with constant K and A(l) = and that the initial measure m 
has a finite second order moment. Then, \/i £ {1, . . . , n}, Vt > 0, 



E[(xi-Yrr]< 



i 



r? 



2+2 



G 



+ E[(A -x) 2 ] + <T 2 t + 2 / / A(F s (x))dx(is 



where X % is solution of (10). 



Proof. Denoting X?(t) = (X 1 , . . . , X t n ), Xf' n (i) = (A t 1,n , . . . , A"' n ) and 
yr(i) = (^,...,n i onek S 

(28) \X?(t) - Y^ n (t)\ 2 = \X?(t) - PX™' n (t)\ 2 

= \X?(t) - PX?(t)\ 2 + \PX?(t) - PX?' n (t)\ 2 

i / n \ 2 n 

(29) <±lj2(Xi-x)) +^(xi-xn 2 . 



n . 

\i=i / i=i 



Since (X t - x) 2 < 3((A - x) 2 + cr 2 ^ 2 + || A'W^t 2 ), the variable X t is square 
integrable. As 

Vx > \(x- x)A{F t (x))\ < P'HqoCI - F t (x))(x + |x|) 

'E(X 2 ) 



< IIA'I 



+ \x\(l-F t (x))) i 
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one has lim x _ >+oc (x — x)A(F t {x)) = 0. Similarly (x — x)A{F t {x)) also van- 
ishes as x — > — oo and J R (x — x)A' (F t (x))pt(x) dx = — J R A(F t (x)) dx. Com- 
puting (Xt — x) 2 by Ito's formula and taking expectations, one deduces that 

E({X t -x) 2 )=E((X -x) 2 ) + a 2 t + 2 f f A(F s (x)) dx ds. 

JO JR 

Moreover, by (3), E,(X t — x) = —A(l)t = 0. One concludes by taking expec- 
tations in (29) then using Theorem 1.5 and exchangeability of the particles. 

□ 

Let us now study the long time behavior of the projected particle system. 

Theorem 2.12. Assume that the function A is uniformly convex on 
[0, 1] with constant a [see (11)] and such that A(l) = 0. Then, the probability 
measure with density 

p n oo(y) = ^e- 2 ^li a -^ 

with respect to the Lebesgue measure dy on A4 n is invariant for the projected 
dynamics (27). Here < y(2) < • • • < y( n ) denotes the increasing reorder- 
ing of the coordinates of y = (yi, . . . , y n ) and Z n = fj^ n e ^ ""W^W &y. 
Moreover, if (Yq 1 ' 11 , . . . , Y n ' n ) admits a symmetric density Pq (y) with respect 
to the Lebesgue measure on M. n , then for all t > 0, (Y t ' n , . . . , Y[ l,n ) admits 
a symmetric density Pt{y) which is such that 

(30) JMAP °° 1 

JMn \P^O J 

2 

where the sequence (X n )n is bounded from below by 12 " .^ • 

In order to deduce long time properties of the nonlinear process from 
long time properties of the projected system, it is not restrictive to assume 
that p^ is symmetric (see Remark 2.15 to get some intuition about this 
hypothesis). But the lack of uniformity in time of the estimation given in 
Proposition 2.11 is a real problem. 

2 2 12 

Remark 2.13. In case n = 2, the process Yt = Y t ' — Y t ' solves the 
stochastic differential equation 

dY t = a(dB 2 - dB\) - sgn(y 4 )(a 2 (2) - a 2 (l)) dt 
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and the density of Y t converges exponentially to a 2( 2 ^~°a(i) e (-(a 2 (2)-a 2 (i))/o- 2 )|i;| 
when the density of Yq is close enough to this limit. As {Y^ 2 ,Y 2 ' 2 



x + 7j(— Zt, Zt), one easily deduces exponential convergence of the density of 
(Y t ^ 2 ,Y 2 ' 2 ) on the straight line M 2 to ^ffl e-^^wX 1 )/' 8 )^) 



The proof of Theorem 2.12 relies on the following Poincare inequality. 

Proposition 2.14. Under the assumptions of Theorem 2.12, the den- 
sity 

p^y) = nll {yi<V2<-<yn} c -(2/a 2 )Y* ^fflw 

Z n 

on Ai n is such that for f :M. n — ► R regular enough, 

f(y)- I f(y)pUy)dy] 2 p n 00 {y)dy 

- Mn \ JMn / 

(ol) 

T 2 



<f / |PV/(y)fe(y)dy 

JAI„ 



where the sequence (A n ) n is bounded from below by 12 " (j2 . 

PROOF of Theorem 2.12. Let us first check the following Green for- 
mula: for / : R n — > R and u : M n — ► R n regular enough, 

(32) / /V • (Pu) (y)dy = - [ PVf- {Pu) (y) dy. 

Let 1 S R n denote the vector with all coordinates equal to 1 . For <p : R — > R 
and u : M. n — > R n ', one has 

(p(y/nz) f V • (Pv) (y + J dz 

JMn \ V H / 

= / <p(x\ + • • • + x n — nx)V ■ (Pv)(x) dx 

<p'(xi + h x n — nx)l ■ (Pv)(x) dx = 0. 



The function ip being arbitrary, one deduces that J M V • (Pv)(y) dy = 0. 
Since V • P{fu) = Vf • {Pu) + /V • (Pu) = PVf • (Pit) + /V • (Pu), (32) 
follows for the choice v = fu. 

By weak uniqueness for (27), when (Y 1,n , . . . , Yq' 71 ) has a symmetric den- 
sity Pq with respect to the Lebesgue measure on J\4 n , the particles Y t,n , 
iG {!,...,«} are exchangeable and for each t > 0, (l^ 1,n , . • • , l^ n ' n ) has a 
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symmetric density p™. By composition with the projection P, one obtains 
an extension of p™ on R n that we still denote by p™. Since Y^i=i a n{i) = 
n(A(l) - A(0)) = 0, setting 



b (y) - Yl 1 {yT(i)<y T (2)<-<yr(n)} 

r£S n 



/ca n (r- 1 (l))\ 
an(r" 1 (2)) 



one has P6 = b and the infinitesimal generator associated with (27) is Lip = 
?j-V • (FVip) - Pb ■ Vi/j. Computing dijj(Y^ n , Y t n ' n ) by Ito's formula and 
taking expectations then using (32), one obtains 



Mr 



^(y)dtPt(y)dy 



Mr 



Lip(y)Pt(y)dy 



Mr 



^(y)V-P(-Vp^ + b P 2)(y)dy. 



Hence the densities solve the Fokker-Planck equation 

a 2 



« = V-P( y Vpr + 6p t n 



Now using (32) and b= — % V ,f°° , one deduces 

JMn \PoC J 

= 2 / PL(y)V-p(^Vp? + bpA(y)dy 

JMn Poo V / / 



(33) 



-a 



-a' 



JMn Pnr, 



Mr 



P'i 

p n J ' 



-(y)p2o{y)dy 



P™(y)dy. 



By symmetry of the function and (31), 

PoG 



a' 



1 




'M n 


Plo 



p^(y)dy = a 2 

> \ 

> A 



Mr 



P^^r(y) 



pZo(y)dy 

2 



JMn \Poo / 

oo / 



El 

Mn \PcO 
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Notice that the computation in (33) is formal and can only be justified 
when pf is a smooth solution of the Fokker-Planck equation. 



Remark 2.15. Let us denote by Y t 



(l),n 



<■■■<¥< 



(n),n 



reordering of (Y t 



l,n 



the increasing 



, Y™ ,n ). According to [9], the reordered system is a 



diffusion process normally reflected at the boundary of the closed convex set 
{y £ M n ■ Vi < V2 < ■ ■ ■ < Vn}- More precisely, 



(34) 



dY 



(i),n 



adpi-a n {i)dt + (^ t 



(7i-7i +l )d\K\ s ,l<i<n 



/ t>0 

with finite variation equal to \K\t, 



Tt")d\K\ t , 

is a continuous process 



7 1 = 7 n+1 = 0, 



t d\K\ t a.e. V2 < i < n,7| > and ~/j{Y t 



(i),n 



where (/3 , . . . , (3 n ) is a Brownian motion such that 



v (i-l),n 
1 t 



If the initial condition (y Q (1),n < ■ • ■ < Y c 



r(n),ni 



) = o, 



admits a density with 



(l),n 



respect to the Lebesgue measure on A4 n , then the law of (Y t 
is the image by increasing reordering of the symmetric law of the solution 



y-(rc.),riN 



> ± t 



l.n 



,Y t ' ) to (27) starting from (Y Q 



l,n 



(l),n 



, y o ' ) with density j>q ob- 
Y^' n ) has the density 



tained by symmetrization of . Therefore (Y t 
Ptiv) =n] -Pt{v)^{yi<-<y n } an< ^ (30) holds with p n replaced by p n . 

In order to prove Proposition 2.14, we take advantage of the specific form 
of the density p^ . Remarking that is the density of the image of a vector 
of independent exponential random variables by a linear transformation, one 
first obtains the following result. 



Lemma 2.16. The Poincare inequality (31) holds with the constant X n 

2 ~ 

greater than multiplied by the smallest eigenvalue X n of the (n — 1) x 
(n — 1) matrix Q n defined by VI < i,j <n — 1, Q™- = b n (i) Lfjb n (j) where 



i(n-i) 

b n (i) = — and L n 

n 



( 2 - 1 









° \ 


-1 2 


-1 










-1 


2 


-1 


... 





... 





-1 


2 -1 





... 







-1 2 


-1 


V ... 






-1 


2 / 
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The last statement in Proposition 2.14 then follows from the next lemma 
which is obtained by interpreting Q n as a finite element rigidity matrix 
associated with the operator —x(l — x)d xx (x(l — x).) acting on functions on 
(0, 1). The Hardy inequality stated in Lemma 2.18 ensures that it is enough 
to bound the smallest eigenvalue of the corresponding mass matrix from 
below. The resort to this one-dimensional Poincare-like inequality in order 
to estimate the constant in the re-dimensional Poincare inequality (31) is 
striking. 

Lemma 2.17. The sequence (A n ) n is bounded from below by 1/(16 x 27). 

Proof of Lemma 2.16. Let / be such that J M f(y)p^ a (y) dy = 0. 
Since the left-hand side in the Poincare inequality (31) only depends on 
the restriction of / to Ai n , one may assume that \/x G K n ,/(x) = f(Px), 
which ensures that for (x\, . . . , x n ) £ W 1 such that x± + ■ • ■ + x n = 0, f(x + 
xi,...,x + x n ) = f(xi,. . . ,x n ) and PVf(x + xi, . . . , x + x n ) = Vf(x 1 ,...,x n ). 
Therefore the Poincare inequality (31) is equivalent to /(/) < f--f(|V/|) 
where 



Us) 



(9 P^o)(-(x2 + --- + x n ),x^)dx^ with X% = (X2,...,X r , 



To integrate the coordinates over independent domains, we make the change 
of variables = Mx^ where 

1 

1 



/ 2 



M 



-1 






V o 



l 







1\ 






-1 1 







One easily checks that for 2 < i < n, z-i + • • • + z. 

deduce that (re — l)z2 + (n — 2)2:3 H h 2z n _i + z n 

fore 



M~ 



and denoting 



1 

n 



/I 
1 
1 

1 

Vi 



2-re 
2 
2 

2 
2 



3 — n 
3 — re 
3 

3 
3 



4 
4 
4 



re 
11 
n 





1/ 

= x 2 H h x n + Xi and 

~-n(x2-\ \-x n ). There- 



1 \ 

■1 

1 



n 



n 



1/ 



n 



n 



N 



n 



n 



2 
n 



1 

n 



M~ 
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one has 



where 



^nJ(M.+ ) n - 1 \M\ 



Pn(i) = -[(*- l)(on(«) H 1" a n {n)) 

n 

- (n + 1 - z)K(l) + • • • + a n (i - 1))] 
= -n^((i - l)/n) > 0. 

Here |M| denotes the determinant of the matrix M; it is equal to n by an 
easy computation. The one-dimensional exponential density with parameter 
c satisfies the Poincare inequality with optimal constant 4/c 2 . Tensorizing 
this inequality (see Chapters 3 and 6 in [1] for further details), one obtains 



Since ^4 is uniformly convex with constant a and -A(O) = A(l) = 0, 

„ / x . ,, s , x to i — 1 / i — 1 \ a , , 

0n(i) = -nA((i - I) n) > -— X 1 = -b n (i - 1 . 

2 n \ n J 2 



„ n n 4 

/ E E ^AVi^i- A/^/i&(-(x 2 + • • • + x n ),x%) dx n 2 . 

jRn 1 k,l=lj=2 Pn(J > 



Therefore 

4 . n n— 1 



^ 4 . a n—± jy 

/(/)<—/ £ Ewv^/^-(-(^ + --- + ^),^)^ 



<-^-I(|V/|) 
crA n 

where A n denotes the inverse of the largest eigenvalue of the symmetric pos- 
itive semidefinite matrix NN* defined by 2Vjj = & & . To prove Proposition 
2.14 with a possibly modified lower bound, it is enough to check that the 
largest eigenvalue is bounded from above uniformly in n. Unfortunately, 
the trace of the matrix can be bounded from below by a positive constant 
multiplied by log(ra). Therefore one has to be more precise. 

Let w be an eigenvector associated with the largest eigenvalue: NN*w = 
j-w. Of course N*w is nonzero and multiplying the previous equality by 

An _ _ _ 

N* , one obtains that N*w is an eigenvector of N*N associated with the 
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eigenvalue By symmetry, ~— is also the largest eigenvalue of N*N. We 
are going to check that the latter matrix is invertible with inverse equal 
to Q n in order to conclude the proof. Because of the definition of N, it is 
enough to check that N*N is invertible with inverse equal to L n . 

By construction of the matrix N, for the equation Nz^ = x where x £ W 1 
to have a solution z%, it is necessary and sufficient that x\ = — (x2 + - ■■ + x n ) 
and then — Mx% ■ 

Now for fixed y G M n_1 , let us find x% 6 M n_1 such that N*x = y where 
x = — (X2 + • — h x n , ^2). This equation writes 

/ / N u N n ... N n \\ 

N 12 



(M 



\ 



N l2 



\iVi n _i iV lft _i 



N 12 



N ln -J 



Xo 



One easily checks that the (n — 1) x (n — 1) matrix in the left-hand side is 
equal to 



/I 






1\ 

1 
1 





Vo 



1 

1/ 



with inverse R ■ 



/I 








Vo 



-1 







1 



1 






1 






\ 






-1 
1 / 



Combining xJ, = Ry with the solution of the previous problem, one obtains 
that the unique solution of the equation N*Nz% = y is Z2 = MRy. One 
concludes by checking that the matrix MR is equal to L n . □ 



Proof of Lemma 2.17. For i e {1, 



0. 



udx) 



i(n — i)(x — (i — l)/n) 



,n — 1}, the functions 

i - 1 i+ 1 

if xe (o,i)\ 

"i - 1 



n 



^nx(l — x) 
i(n — + l)/n — x) 



nx(l — x) 



if x e 

if x G 



n n 
i i + l 

n ' n 



are such that 

Vi,j G {!,..., n- 1} 



Qij 



(x(l — x)ttj(x))'(x(l — x)uj(x))' dx. 



By the Hardy inequality stated in Lemma 2.18 below, the smallest eigenvalue 
of the matrix Q n is greater than the smallest eigenvalue of the (n — 1) x (n — 
1) tridiagonal matrix Rfj = Jq 1 Ui(x)uj{x) dx divided by 16. 
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For 7 £ {1, ... ,n - 2}, let rf = «i(i*t - u i+1 )(x) dx and 

Z 1 2 / x A _ in -If f 1 1 _ n-l 

J(n-l)/n n J(n-l)/n x n 

Using the change of variables y = 1 — x, one easily checks that 

Vz £ {1, . . . , n - 1} Rl - Rl^ - Rl +l = r? + rJU 
where by convention i?" = = 0. We are going to prove that 

Vn>3 Vi£ {2,..., n -3} r?>^, 
and that r" and 7"^_2 are nonnegative. For y £ IR™" 1 , one deduces that 

n— 1 n— 2 

7/*iA/ = 53 i^ y 2 + 2 53 

i=l i=l 

n-1 n-2 , ,2 

i=l i=l 

and the conclusion follows. 

Let us first suppose that i < |_§J — 1, which ensures that the function 
f(x) = x 2 (l — x) 2 is increasing on [i/n, (i + l)/n\. Let g(x) = Ui(m — Ui+i)(x). 
One easily checks that 

( i+1 )/ n . . , i 2 (n-i) 2 (\ + 

g(x) dx = 

i/n n 4 \3 Qi(n — i) 

0, if 7 = 1, 

> { i 2 (n-i) 2 



12n 4 



if i > 2. 



Since there is some X\ £ [i/n, (i + l)/n] such that the function g{x) is non- 
negative on [i/n,Xi] then nonpositive on [xi, (i + l)/n] , and / is positive and 
increasing, one deduces that for all x £ [i/n, (i + l)/ra], J^ n jM dy >0. This 
ensures that Vx £ [i/n, (7 + l)/n] 

f/w r fri^) =/'(*) r f^^+^)> 9 (x). 

axV Ji n fiy) J Ji n fiy) 



Therefore 



Ai+l)/n a ( v ) i r(i+l)/n 

r?= jr\ dy > 777— -ttt-t / 9 (y)dy 

Ji/n f{y) flit + l)/n) Ji/„ 

0, if 7 = 1, 
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Let us now suppose that i > \ so that the function / is decreasing on 
[i/n,(i + l)/n]. We deduce that 

1 

r?>—r- (fu 2 )(x)dx 



f(i/n) Ji/ 



1 /•(*+!)/" 



/((i + l)/n) 7 V 
1 z(n — i) 



3 6(* + l)(n-i-l) 

and the left-hand side is greater than 1/12 for i < n — 3 and nonnegative for 
i = n — 2. 

We still have to deal with the case n odd and i = (n — l)/2. Then, / is 
not monotonic on I n = [i/n, (i + l)/n] = [1/2 — l/2n, 1/2 + l/2n]. But by 
symmetry, 

1 \2 /•l/2+l/2n /l /o i 1 /% _ ^Vl _ Oo."\ 

C?X 



(n — 


l) 2 (n 


+ 1) 2 




16n 




(n — 


l) 2 (n 


+ 1) 2 




32?i 




(n — 


l) 2 (n 


+ 1) 2 


2n 



(n ' 1)/2 16n A/2-l/2n " X 2 (l-x) 2 

■ dx 



l/2-l/2n X 2 (l-X^ 



2 



2+l/2n ( 2_i\2 

> ^ v " ' ^ / (l-2x) 2 dx^ (n 



l/2-l/2n 



6ra 4 



which completes the proof. □ 



Lemma 2.18. For a// «£ L 2 (0,1) such that the distribution derivative 
(x(l -x)u(x))' belongs to L 2 (0, 1), 

u 2 (x)dx< 16 / ((x(l — x)u(x)) ) dx. 
Jo 

Proof. For v a C°° function with compact support on (0,1), by the 
integration by parts formula, 

1/2 V 2M f l/2 v 2, x ) j rl/2 vv '( x \ 

o x 2 (l — x) 5 



dx<4 —±-Ldx = 8[ y —^dx-v z {l/2) 



r-l/2 ^(V) \ 1/2 / rl/2 n \ 1/2 

< 8 (/ ^dx) (/ (v'(x)) 2 dx^ 



x* 

Dealing with the integral on (1/2, 1) in a symmetric way, one deduces 

(35) ^ J' (X) N9 dx < 16 l\v'{x)) 2 dx. 

Jo x z (l-x) z Jo 
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Now approximating v € Hq(0, 1) by a sequence of C°° functions with com- 
pact support converging in the H 1 norm and almost everywhere, one deduces 
with the Fatou lemma that the inequality still holds for v € Hq . 

For u satisfying the hypotheses in the lemma, v(x) = x(l — x)u(x) be- 
longs to i/ 1 (0, 1). According to Theorem VIII. 2, page 122 of [4], v admits 
a representative continuous on [0, 1] still denoted by v. Moreover, since 
u(x) = j^-ffj.) belongs to L 2 (0, 1), necessarily, v(0) = v(l) = 0. By Theorem 

VIII. 11, page 133 of [4], v belongs to Hq(0,1) and the conclusion follows 
from (35). □ 
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